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Abstract. In this small review we present the actual state the knowledge about weighting 
black holes. Black holes can be found in stellar binary systems in our Galaxy and in other 
nearby galaxies, in globular clusters, which we can see in our and nearby galaxies, and in 
centres of all well-developed galaxies. Range of values of their masses is wide and cover 
about ten orders of magnitude (not taking into account the hypothetic primordial black 
holes). Establishing the presence of black holes, and in particular the measurement of their 
mass is one on the key issues for many branches of astronomy, from stellar evolution to 
cosmology. 
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1. Introduction 

Astrophysical black holes (BH) are customar- 
ily divided into three classes: stellar mass black 
holes, with masses of a few up to twenty-thirty 
solar masses, intermediate mass black holes 
(IMBH) with masses of a few hundreds to a 
few thousands of solar masses, and massive (or 
super-massive) black holes residing in galactic 
dynamical centres, with masses from 10 5 solar 
masses up. 

Determination of values of black hole 
masses is important for the three main rea- 
sons: (i) to distinguish small mass black holes 
from neutron stars (ii) to analyse the luminos- 
ity states of accreting black holes through the 
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luminosity to the Eddington luminosity ratio 
(ii) to study the growth of central black holes 
in galaxies. All three aspects are important, the 
last one is a key issue in cosmology, in studies 
of the galaxy evolution and AGN feedback. 

The methods used for mass determination 
can be broadly divided into the following three 
classes: 

- dynamical methods 

- spectra fitting methods 

- scaling methods 

although such a division is not necessarily 
unique. For example, the reverberation method 
is rather an example of a dynamical method 
while the Broad Line Region (BLR) radius 
vs. luminosity relation combined with a line 
width fitting is somewhere in between the di- 
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rect dynamical method and a simple scaling 
law. In this review we will include it as a scal- 
ing method. 

There are several recent reviews covering 
the topic of black hole mass determination 
(e.g.lCasaresl2007t Ziolkowski 1120081: IWandeTI 
l2008tlVestergaardll2009l) . Here we will not go 
into the details of the various methods and its 
caveats but instead we will provide the collec- 
tion of laws with short description of applica- 
bility. The possibilities are broad so black hole 
mass estimate is possible nearly for every ob- 
ject. Different methods use observational data 
from various energy bands from radio through 
NIR and optical to X-ray domain. Therefore, 
the black hole mass measurement is one of the 
nice examples of the multifrequency approach. 

2. Black hole mass determination 
methods 

2.1. Dynamical methods 

Dynamical methods are mostly based on the 
Keplerian motion of a test particle around the 
black hole. The basis idea is simple: if for a 
test particle we can measure the velocity v at 
the circular orbit and the radius R of the orbit 
we obtain the value of the black hole mass: 



M = 



v z R 



(1) 



The application is complicated by several fac- 
tors. First, we have to find a suitable test parti- 
cle. Next, there is a problem of the inclination 
of the orbit. Even for a circular orbit the mea- 
surement of one of the velocity components 
does not allow for a unique determination of 
the orbital speed. 

We have a choice of the following test par- 
ticles and measurements: 

- Companion star (in binary systems) - mass 
function 

- Companion black hole (e.g. OJ 287) - or- 
bital period 

- Nearby star (e.g. S02 star in Sgr A*) - ap- 
parent motion, radial velocity 

- Nearby stars - stellar rotational velocity, 
stellar velocity dispersion and/or stellar lu- 
minosity profile 



- Broad Line Region clouds - line width and 
BLR radius from reverberation 

- H20 masers - velocity and radius from 
mapping 

Of those methods the last one - water maser 
method - is the most accurate since we ob- 
serve water maser emission in an AGN only 
when observing almost in the orbital plane of 
the emission. Additionally, since the emission 
of several blobs located at various distances 
from the black hole is measured we have a 
test whether the motion is indeed Keplerian. 
The emission comes from the geometrically 
thin disk, from distances of order of less than 
a fraction of a parsec from the central black 
hole. Such measurement w a s done for NGC 
4258 dMiyoshi et al. 1 1 19951: iHerrnstein et al. I 
120051) . and the most recent measurements give 
the black hole mass value of (3.82 + 0.01) x 
1O 7 M water maser (see ISiopis et al. I I2009L 
and the references therein). Similar measure- 
ments were done for a few other sources, but 
the water maser emission is rare (detected in 
about 10 percent of the nearby Seyfert 1 .8 - 
Seyfert 2 galaxies, Kon dratko et al. 1120061) . as 
it requires special inclination and relatively low 
absorption as for a Seyfert 2 galaxy (highly 
inclined AGN are extreme type 2 AGN), the 
data quality is usually not good enough to al- 
low for mass measurement (emission is unre- 
solved), and in some cases like in NGC 1068 
the implied motion is clearly non-Keplerian so 
the mass measurement cannot be performed re- 
liably. The method does not apply to Seyfert 
1 galaxies or lower mass black holes, and the 
majority of Seyfert 2 galaxies also do not show 
maser action so the method has very limited 
applicability. 

In case of stellar mass black holes, since 
they are generally in binary systems we have a 
secondary playing the role of the test particle. 
The companion star, however, is itself massive 
so the measurement of the companion semi- 
amplitude of the radial velocities, K c , and the 
period, P, of the orbital motion results only in 
determination of the mass function 



f(M) 



(M sin i c y _ KlP 
(AT + M c f ~ 2^G' 



(2) 
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which gives only the lower limit for the black 
hole mass M. The further determination of 
the black hole mass requires both an estimate 
of the inclination of the orbit, z' c , as well as 
the estimate of the mass of the companion 
star, M c . Best constraints are achieved if mul- 
tiwavelength measurements are used (e.g. to 
confirm the binary period) and evolutionary 
constraints applied. New determinations of the 
black h oles were recently performed for exam- 
ple bv lOrosz et alTI (l2009h (LMC X-l, M = 
10.91 + 1.41M ). 

Overall, black hole masses in binary sys- 
tems of the Milky Way were determined for 
24 objects (confirmed BH). The accurate val- 
ues (err or less then 10%) were obtained for 7 
objects dCasare s 2007; Ziolkowski 2008). 

Similar analysis of the binary motion was 
also used in the case of two massive black 
holes. In OJ 287 the observed periodic out- 
bursts every 12 years are interpreted as a 
result of a smaller massive black hole (~ 
1.3 x 10 8 M o crossing the accretion disk of a 
mor e massive (~ 1.8 x 10 10 M o ) black hole 
(see IValtonen et al. 1 120091 and the references 
therein). However, this object is exceptional, 
with the optical data covering more than 100 
years. The search for other massive binary 
black hole cont inues but the list of ca ndidates 
is short (see e.g. lXu & Komossa 120091) and the 
available data do not allow for black hole mass 
measurement at the basis of binarity itself. 

Single star as a test particle can also be 
used in the case of the massive black hole 
mass measurement for Sgr A*. In this case the 
star is indeed a test particle, but stellar orbits 
are elliptical and inclined. The peculiar mo- 
tion of these young bright stars can be fol- 
lowed using the active optics in the IR band. 
One of the stars (S02) already completed the 
whole orbit and the most recent determination 
of the Sgr A* mas s based on that result is 
(4.1±0.6)xl0 6 M o (Idhez et al. 120081) . Similar, 
and even more accurate value is obtained if all 
28 mon itored stars are used ( (4.3 1 + 0.06) x 
10 6 Mm. iGillessen et al. I l2009h . The provided 
uncertainty of the measurement is addition- 
ally enhanced to 0.36 x 1O 6 M due to system- 
atic error in the estimate of the distance to the 
Galactic centre. 



The method based on a single star cannot 
be used now for other, even nearby galaxies or 
for centres of the globular clusters (due to the 
crowded fields). 

Thus the next approach is to use the mea- 
surement of the rotational velocity or/and ve- 
locity dispersion of the stellar distribution to 
infer the presence of the central black hole. The 
method is not a simple one since it requires 
modelling of the stellar motion on the top of 
the measurement of the stellar velocity from 
the width of absorption lines. If the measured 
stars are not close enough so that the gravita- 
tional potential contains the contribution from 
the stars themselves the problem is particularly 
difficult. 

The black hole mass measurement for 
nearby non-active black holes based on 
stellar dispersion was performed for several 
objects observed by HST or ground-based 
telescopes, preferably with active optics 
(iMagorrian et al. II 19981: iGebhardt et al. Il2000t 
TloOOt iKrajnovic et al. I 
shown for example by 



Ferrarese & Merritt 
2009) 



As w as 

iGiiltekin et al~l d2009al) for NGC 3607, the 
fit with a central black hole of mass of 
(1.2 + 0.4) x 10 8 M Q is better by X 2 = 10.6 
than the fit without a central black hole. 
Application of the method to globular clus- 
ters in our Galaxy is even more difficult 
since the size of the region dominated by 
the central black hole has smaller angular 
size. Minor changes about the assumption 
of the stellar dynamics change the result 
drastically. For example, the fit to the globular 
cluster u> Centauri gives no black hole for 
core assumption for the stellar distribution 
but for the cusp profile assumption it gives 
the black hole mass of (8.7 ± 2.9) x 1O 3 M 
(Ivan derMarel& Anderson I [2009). Another 
ex ample of the potentia l problems was shown 
bv lFregeau et al. I d2009l) : the inclusion of natal 
kicks of white dwarfs in globular cluster NGC 
6397 removes the need for the IMBH at the 
cluster centre. 

The determination of the black hole mass 
directly from the stellar dispersion velocity 
profile should not be confused with a global 
measurement of bulge velocity dispersion and 
black hole mass determination based on the 



Czerny & Nikolajuk: Mass of black holes 



285 



scaling between the two since the motion of 
the stars in bulge is not determined dynami- 
cally by the black hole mass; this method will 
be described in Sect. 12.31 

Dynamical methods can be extended to dis- 
tant objects for active galactic nuclei since in 
this case the gas of the Broad Line Region 
(BLR) can play the role of the test particle. 
The BLR is never resolved but the emission 
line profile provides the information on the ve- 
locity and the orbital radius can be measured 
from the time delays between the variations of 
the continuum and the resp onse of the line, i.e. 
from reverberation (see e.g. lPeterson & Home I 
|2004| for the general description). 

The line profile is converted to the velocity 
either through measurement of the full-width- 
at-half maximum (FWHM), or through line 
dispersion (seco nd moment of the line profile), 
v a . As argued by lDenney et al. 1 (120091) . the first 
method is better for low signal to noise spectra 
while the second method is more accurate for 
higher signal to noise spectra since it is then 
less sensitive to details of the subtraction of the 
narrow component of the line profile. 

The black hole mass is then determined 
from the relation: 



Mi 



BH 



fl 



D,.2 

FWHM 

G 



or from the relation: 

Rvl 



(3) 



(4) 



The problem of the method is the geom- 
etry of the BLR. If we assu me the spheri- 
cal geometry, as introduced by lNetzerl (119901) . 
the coefficient / = f in the formula above 
is deter mined uniquely. This approach wa s 
used bv IWandel. Peterson. & Malkanl (1 19991) : 
iKaspi et al. I d2000t) . 

How ever, the BLR is likely to be flattened 
(see e.g. iMcLure & Dunlopl bOQl) so the co- 
efficient was later determined at the basis of 
the comparison between two mass determina- 
tion methods (reverberati on and stellar disper- 
sion, lOnken etal. 1 120041) and this coefficient 
was adopted in Eqs. (01 and (@). It is by a fac- 
tor ~ 1.8 higher t han in spherically symmetric 
case. In addition. [Collin et aTl d2006l) reported 



that the value of the coefficient / is dependent 
on type of AGN (for example f\ in Eq. <£3j for 
Seyfert 1 galaxies with broad emission lines 
is less than in case of Narrow Line Seyfert 1 
(NLS1) galaxies). 

Reverberation is time consuming but nev- 
erth eless it was performed for over 40 objects 
("see iPeterson et al. 1 120041 iBentz et alJ 12009 eL 
Bent z et al. 12009 cj, and the references therein). 
It includes several Seyfert 1 galaxies, about 15 
quasars and 7 NLS1 galaxies. 

The problem in this method is the poten- 
tial dependence of the scaling factor on the 
inclination angle. The iss ue was discussed in 
a number of papers ( e.g. ICollin et al. I l2006t 
iNikoiajuk et al. 1120061) . and the systematic er- 
ror even by a factor of 3 can arise du e to this 
effect in some cases (see discussion by lKrolikl 
J2001K 

Additionally, iMarconi etaTI d2009l) sug- 
gested that the effect of radiation pressure 
should be included in determination of the 
black hole mass based on dyn amical methods 
(i.e. virial theorem). However. lOnken I d2009h 
argued against the importance of this effect. 

2.2. Spectra fitting methods 

The spectra fitting methods are based on the as- 
sumption that we have a reliable model of the 
emission coming from the vicinity of a black 
hole. Such a model usually depends on the 
black hole mass, but also on some other param- 
eters, mostly accretion rate, but also the spin of 
the black hole and the viewing angle of the nu- 
cleus. Therefore, using those methods is rather 
difficult and the results should be taken with 
much care. 

In high accretion rate sources, like Galactic 
black holes in their soft states or Narrow Line 
Seyfert 1 galaxies and quasars, the emission is 
mostly dominated by the emission of an accre- 
tion disk seen in X-rays (GBH) or optical/UV 
(AGN). In exceptional AGN the peak emis- 
sion is also seen in X-rays (RE J1034+396, 
ISoria & Puchnarewicz Il2002l) . 

If we have a single measurement of the 
monochromatic flux at vo (in the power-law re- 
gion of the big blue bump) and the position of 
the maximum of the disk component spectrum 
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on vF Y diagram, v max , then assuming a non- 
rotating black hole, disk inclined at 60 deg and 
extending to the innermost stable circular orbit 
(ISCO) and a local black body without any cor- 
rection connected with the colour temperature 
we can have a black hole mass estimate: 

log M BH = - log L Vo - - log y + 
2 O 

4 

-J log Vmax + 16.515 (5) 



dTripp et al. I [[994). This can be applied to 
bright quasars if only broad band low qual- 
ity spectrum is available. In Seyfert galaxies 
the spectral slopes are much softer than ex- 
pected from the Shakura-Sunyaev disk emis- 
sion (i.e. the effects of X-ray irradiation, ex- 
tinction and outflow are important) so the for- 
mula is inadequate. For higher quality broad 
band spectr a more advanced models can be fit- 
ted (see_e^_ Davis. Woo. & Blaes1 f2007 based 
on iHubenv et al. 112000 models; irradiation ef- 



fect in Loska, Czernv, & Szczerba 2004 and 
ICackett et al. 112007b ICzerny & Janiuk II2007I) . 

For galactic sources, mor e advanced disk 
spect ra model are used (e.g. Sadowsk i et al. I 
2009) and with the knowledge of the black hole 
mass from dynamical method (binary motion) 
they are used instead to infer th e black hole 
spin (e.g. lMcClintock et al. f2006l) . 

However, in order to rely on such mass 
determination we must be sure that the mea- 
sured emission does come from a fairly stan- 
dard accretion disk. In some sources, like 
ULX, this is highly uncertain. For example, 
thermal disk + pow er law fits to X-ray spec- 
tra im ply IMBH (e.g. lMiller. Fabian. & M iller 
2004) while other decomposition of the spec- 
tra implies the ultral uminous high states of 
100 Mo BH (e.g. iDone & Kubotal l2006t 
iRoberts 1 120071) . Recent study of the sample 
of 94 Chandra sources dBerghea et al. 1 120081) 
shows that the cool disk spectral component 
is unrelated to luminosity which implies that 
indeed we do not understand the spectra! 
Arguments against low masses in ULX sources 
(i.e. against strong beaming) are rather indi- 
rect (see the discus sion of the sourc es ULX- 
1 ESO 243-49 by (IFarell et al. 1120091) for a re- 
cent example of such analysis). Therefore, con- 



clusions on the presence or absence of IMBH 
based on spectral analysis cannot be done. 

Accreting black holes also show the pres- 
ence of the X-ray emission and the reflection 
component coming from X-ray reprocessing 
by an accretion disk. The iron line formed in 
this process can be also used for measurement 
of the black hole parameters. The line profile 
is mostly sensitive to the black hole spin but if 
X-ray time variability is included then iron line 
reverberation also should allow (in principle, 
when next generation large area X-ray tele- 
scope will operate) for black h ole mass mea- 
surement dDovciak et al. 112004 . 

2.3. Scaling methods 

Scaling methods are invaluable for black hole 
mass determinations in large AGN surveys as 
well as when a simple black hole mass es- 
timate is needed for a specific object. There 
are numerous methods based on various ob- 
servables. The accuracy of the methods is dif- 
ficult to assess for a single object so two or 
more methods should be used whenever pos- 
sible. These methods are based on scaling of 
measured properties with the mass found ear- 
lier by independent mass and the discuss quan- 
tity measurement. 

Most broadly used nowadays is the method 
of the black hole mass measurement from a 
single spectrum in IR, optical or UV band. 
The method originally comes from reverber- 
ation approach but is based on empirical 
connection established between the BLR ra- 
dius and th e source monochromatic lum i nosity 
(see e.g. IWandel. Peterson. & Malkan 1 1 19991: 
iKaspi et al. Il2000h . If the starlight contamina- 
tion is likely to be unimportant, or it is ac- 
counted for, the best relation to use is 



log 



R 



BLR 



1 It day 



-21.3 + 



+(0.519!°;°|!)log(,lL /l (5100A)) 



(6) 

(based on HST spectra. iBentz et al.ll2009al) and 
consequently the formula for the black hole 
mass determination based on H/3 line reads: 



log M BH = log 



/FWHM(Hy6)\ 2 
I 1000 kms" 1 / 
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( AL A (5100A) V 
\ lO^ergs- 1 / 



0.519 



+ log/i +6.81 



(7) 



where f\ is the factor from Eq. (O and Mbh is 
in units of M e . If the starlight contamination 
is likely, the older version of the phenomeno- 
logical fit may be more accurate in practice. 
iKaspi et alTI (|2005) give the relation 



10 It days 



2.39 



AL A (5100A) 
lO^ergs- 1 t 



,0.67 ±0.05 



(8) 



and from this it follows that 

/FWHM(H/J)\ 2 



log M BH = log 



/ /tL A (5100A) V 
X l lO^ergs- 1 j 



\ 1000 kms- 1 / 

0.67" 

+ log/! + 6.66 



(9) 



In their A GN analysis 

IVestergaard & Peterson ] (|2006) use two 
formulae: 



log M BH = log 

/ ^(5100A) V 
A lO^ergs- 1 / 



/FWHM(H/3)\ 2 
I 1000 kms- 1 / X 



+ (6.91+0.02) . 



(10) 



or 



log M BH = log 



UO^ergs- 1 / 



/FWHM(H/3)\ 2 
I 1000 kms" 1 / 



+ (6.67 + 0.03) . 



(ID 



If H/3 is not available, other lines can be 
used. For example, the recent formula based on 
the Mg II line reads: 



log M BH = log 



/FWHM(Mgn)> 



\ 1000 kms- 1 



I AL A \» 
I lO^ergs- 1 ) 



+ zp(A) 



(12) 



where zp(A) = 6.72, 6.79, 6.86 and 6.96 for 
^1350A, /121 00A, ^3000A, and /1 5100A , 
respectively (IVeste rgaard & Osmer ] l2009h . 



There is also a version of this formula based 
on the line flux instead of the continuum 
mono chromatic flux (see e.g. Ko ng et al. I 
120061) . 

If even Mg n is outside the measured spec- 
trum, there is a possibility to use C IV line 
(IVestergaard & Peterso"n"1l2006l) : 



log M BH = log 



/FWHM(C iv) \ 
I 1000 kms- 1 ) 



q^(1350A) V 
v lO^ergs- 1 j 



0.53 



+ (6.66 + 0.01), (13) 



This line is perhaps less reliable as belong- 
ing to High Ionisation Lines (HIL) and being 
more likely influenced by the radiation pres- 
sure that Low Ionisation Lines (LIL) (like H/3 
and Mgll). 

Those methods work for type 1 AGN, 
i.e. Seyfert 1 galaxies and quasars, in Seyfert 
2 galaxies we do not see the BLR region. 
However, in those sources we can use scaling 
based on the bulge properties which generally 
seem to apply for all AGNs well as non-active 
galaxies. 

The best relation, with very low disper- 
sion, was found between the stellar disper- 
sion in the bulge and t he black hole mass 
dGebhardt et al. I 120001: iFerrarese & Merritt I 



2000). Gultekinelaljj2009b) based on obser 
vations of over 49 observations of galaxies give 
formula: 



log M B h = 8.12 + log 



1 200 kms- 1 / 



4.24±0.41 



(14) 



k 200 kms- 1 , 
which implies the surprising relation be- 
tween the black hole mass and the bulge 
mass dKormendv & R ichstone 1 1 9951: 
iMagorrian et al. I 1 1998h . If thestellar dis- 
persion in the bulge is not measured we 
can use first the relation between the bulge 
luminosity and the bulge mass 



l0g (^_, 11 + U8l0g (^) a „ 

Then we can combine it with the black hole 
mass vs. bulge mass relatio n in order to receive 
M BH (iHiiring & Rix 1120041) : 

Mbulg- ^ 1-12*0.06 



Mi 



BH 



S / A/bulge V- lz±u - uo 
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If the bulge luminosity is measured in B band 
then we have to recalculate it to V band 
through empirical r elation B - V = 0.8 (e.g. 
iBian & Zhao l2003 | ) or directly use the form ula 
proposed by lKormendy & Gebhardtl Q2001 ) 

M » = °- 78 * loS fef V °- <"> 

For AGN the most recen t scaling relation 
found by lBentzetaD (l2009b) is: 

(L V bulge \°- 80±009 
log M BH = 7.98 + log (j^^j (18) 

For highly inclined radio galaxies the fol- 
lowing relation between host galaxy absolute 
magnitude at R-band, M R and M B h can be ap- 
plied 

log M BH = -0.5M R - 2.74, (19) 

dMcLure & Dunlop 120021) where M R is the ab- 
solute R-Cousins bulge luminosity. 

For BL Lac objects the similar relation, 
but with a d i fferen t shift, was proposed by 
iBettoniet all (120031) : 



log M BH = -0.50M R -3.00. 



(20) 



More complex method of deriving the black 
hole mass from the Two-Micron All-Sky 
Su rvey K-band bulge lumi nosity was presented 
by lVasudevan et al. I (120091) . 

Finally, 278 radio-loud AGNs (includ- 
ing 146 BL Lac so urces) were observed by 
IZhou&Cao I (120091) . A significant correlation 
was found between the Lorentz factor of the 
jet, 7 m in, and black hole mass which opens a 
possibility to estimate the last quantity: 



log M BH = 3.26 log y rain + 5.81 . 



(21) 



In case of AGN we have also a possibility 
to use the profile of optical [O in] line (from 
NLR) as a proxy for the stellar dispersion ve- 
locity as advocated by several authors and use 
the Eq. d!41 l. The accuracy is improved with the 
use of a relation 

o% = FWHM([0 m)/2.35 (22) 

dGaskelll [2009b . This possibility is particu- 
larly interesting for radio galaxies where ob- 
scuration of the central region is consider- 
able due to the host galaxy when the source 



is seen in the host galaxy plane, as it is in 
CSO sources. When mid-IR spectrum is avail- 
able, ot her NLR lines can be used, like [NeV], 
[O iv] (iDasyra et al. 112008b . 

Finally, also X-ray data provide us with the 
black hole mass measurement possibility. One 
of the advantages of those methods is the ex- 
pected independence on the inclination effects 
since the X-ray emission is rather isotropic. 

X-ray lightcurves of accreting black holes 
are dominated by aperiodic red noise and 
power spectral densities (PSDs) are well rep- 
resented by several Lorentzians. Such anal- 
ysis is mainly done using high quality X- 
ray data for X-ray Bina ri es in our Galaxy 
(see e.g. iNowak I l2000fc iPottschmidt et al. I 
I2003L for Cyg X-l, GX 339-4). In case of 
low quality X-ray data, generally in AGNs 
and ULXs, PSD of those sources are rep- 
resented by_a__gow^r_Jaw^_wjA_o2ieor_jwo 
breaks dUttlev. McHardv & Papadakis 1 120021: 
iHeil. Vauahan & Roberts 2009). However, re- 
cently |McHjirdy_etay (120071) reported finding 
2 Lorentzian components in PSD of nearby 
NLSl galaxy Ark 564. 

The overall shape of the power spectrum 
depends on the luminosity state (the effect well 
studied in galactic black holes) and shifts lin- 
early with the black hole mass. High frequency 
break was suggested t o scale with mass in 
AGN (Papadakis 2004) as 



1-6 



7 1.7x10" 
M BH = 10 7 M , 

/hfb 



(23) 



wh ere fhfh is the high f requency break in 
Hz. iMcHardv et al. I d2006l) suggest a formula 
which includes also the luminosity term: 

log/ hfb [day] = -2.171og( I ^) + 

+a901og (io^Fr) + 2 - 42 - (24) 

However, the form and strength of the lu- 
minosity dependence is an open issue and 
likely depends of the AGN type (Seyfert 1 
vs. NLSl vs. LLAGN). In galactic sources the 
low frequency part of the X-ray power spec- 
trum clearly shows a dependence on the hard 



Czerny & Nikolajuk: Mass of black holes 



289 



or soft state of the source. No strong variabil- 
ity of the high frequency part of the X-ray 
power spectrum is s een for galactic source s 
in their hard states dGierlinski et al. 1 l2008b). 
Similar disconnection with luminosity and de- 
pendence only on Mbh is see n in Sgr A* 
variab ility observed in Near IR dMeyer et al. I 
120091) : the break in the power spectrum of the 
NIR lightcurve fits expectations from AGN de- 
spite the huge difference in the bolometric lu- 
minosity between AGN and Sgr A*. 

An interesting version of this scal- 
ing emerged from the approach of 
lHavashida et al. I (Qj>98) and uses only the 
high frequency part of the power spectrum, 
or, in practice, the X-ray excess variance 
measured in the timescales shorter than the 
timescale corresponding to 1 //hfb in previous 
formulae. The relation 



Mi 



BH 



1.92- 



T-6t 



-Mr, 



(25) 



dNikolaiuketal. 1120061 [2009) is based on the 
scaling with Cyg X-l and uses 20 M for 
Cyg X-l mass; study for a number of galac- 
tic sour ces gives somewhat lo wer coefficient 
of 1.24 dGierlinski et al. Il2008bl) . Here T is the 
duration of the lightsurve and 5t is the time bin, 
both in seconds, and the X-ray excess variance, 
cr;r ms is dimensionless, normalized by the aver- 
age flux. 

In addition to broad band power spectra, 
quasi-periodic oscillations (QPO) are seen oc- 
casionally in X-ray Binaries and ULXs. Their 
use for black hole mass measurement is com- 
plicated by the fact that there are several types 
of QPO, and firm identification of the QPO 
type is needed to infer the black hole mass. 
QPO speci fic example of the application can 
be found in Shaposhnikov & Titarchuk 1 d2007l) 
who derived the black hole mass in Cyg X-l 
of 8.7 + O.8M from the black hole mass in 
GRO J1655-40 at the basis of their QPO scal- 
ing with mass and X-ray spectral slope. There 
were several claims about the presence of QPO 
also in AGN, and recently iGierliriski et al. I 
(2008a) reported about discovery of an un- 
mistakable QPO in NLS 1 RE J1034+396. 
Middl eton & Done I d2009l) have estimated the 



mass of a black hole in this object (2 - 3 x 
1O 6 M ). 

However, caution in the case of X-ray 
methods is also needed since some weak in- 
clination effect can be expected due to the 
general relativity causing in general slight 
enhancement in variability for highly in- 
c line d objects seen in theoretical studies 
(e.g. Abramowi cz & Baolll994t ICzerny et al. I 



2004). 



3. Consistency checks 

Since all the methods may contain systematic 
errors the best test of the accuracy of the black 
hole mass determinations is to use two or more 
independent methods. This is particularly im- 
portant for all the scaling approaches. In the 
case of AGN, also black hole mass determi- 
nations done at various epochs are interesting 
since they rely on the spectra which vary in 
time. Here we list a number of such compar- 
isons. 

A test of the reliability of a single epoch 
mass measurement based on the BLR radius- 
luminosity scaling (see Sect. 12.3b was per- 
formed using multi epoch measurements for 
two galaxies (NGC 5548 and PG 1229+204, 
iDenney et al. I l2009h . The dispersion of the 
mass measurement due to the source variabil- 
ity was less than 0. 1 dex for high signal to noise 
spectra (above 20). 

The use of vari ous optical/UV lin es is 
nicely illustrated in IVestergaard I d2009l) (see 
her Fig. 3). There are visible discontinu- 
ities in the mass distribution with redshift 
due to the forced change from H/3 to Mgn 
and finally Civ, but less th an a factor 2. 
iRisaliti. Young & Elvis I d2009l) confirm previ- 
ously suggestions (see e.g. IShen et al. Il2008h 
that C iv line is not a reliable indicator of BH 
masses. Much better is to use Mg II. However, 
BH masses determined from this ion has sys- 
tematic error so the Mg n mass can be cor- 
rected using the relation log [M B}1 (H/3)] = 1.8x 
log[M BH (Mgii)]-6.8. 

The comparison of two methods was done 
for a water maser source NGC 4258, with the 
black hole mass of (3.82 + 0.01) x 1O 7 M 
from water maser (IHerrnstein et al. 1 120051) and 
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(3.3 + 0.2) x 10 7 M o from stellar dispersion 
(Siopi set al. 112009b . Stellar dispersion against 
reverberation was also tested for NGC 4151, 
stellar dispersion g ives (4 - 5) x 1O 7 M 
dOnken et al. 1 12007|) and reverberati on gives 
4.57^7, x 10 7 Mr, feentz et al. 112006b 

Based on such tests, IVestergaardl |2009) 
estimates that the systematic error in secondary 
estimates (i.e. scaling laws) is likely to be from 
0.3 dex for stellar dispersion in the bulge to 0.7 
dex when using bulge luminosity or NLR line 
width. 

Interesting tests of the effects of inclina- 
tion in AGN were performed by comparison of 
the reverb eration method and the X-ray excess 
variance (Nikolaiu k et al. 2006). The change 
of the inclination from ~ 15 deg to ~ 60 deg is 
likely connected with the change in the factor 
/ in Eq. [3]by a factor 4 although errors in this 
analysis are large. 

4. Summary of the results 

4.1. Galactic black holes 



The typical accuracy in the determination of 
the black hole masses is about a factor 2 or 
better. The best known black hole mass (mi- 
croquasar GRO J 1655-40, M R h = 6.3 ± 0.5, 
iGreene. Bail yn&Oroszl EOQl) was achieved 
due to the dynamical method (accurate period 
determination, inclination from modelling of 
the lightcurve due to the companion ellipticity, 
mass ratio also from modelling the ellipticity). 
Overall, there are over 40 (24 confirmed BH) 
determinations of the black hole mass in such 
systems in our Galaxy and a few in the bina- 
ries in nearby galaxies like in LMC and M33 
(e.g. eclipsing binary M33 X-7, with black ho le 
mass of 15.65 ± 1 .45M m : lOrosz et al. 112007b . 

4.2. Intermediate mass black holes 

This is the hottest subject in the Art of black 
hole mass measurement. First, the exact range 
is not specified, so it is not clear whether rela- 
tively small but otherwise typical AGN (i.e. lo- 
cated at the centres of their host galaxies) black 
holes from the range ~ 10 s - 10 6 M o (s ee sam- 
ple of iDesroches. Greene & Ho 11200 9) belong 



to this class. T heir existe nce does not cause 
much doubt (e.g. lHo 12008b . The on-going stel- 
lar dispersion measurements in hu ndreds of 
near-b y galaxies (see for example iHo et al. 1 
2009) may bring more of them although mea- 
suring black holes with masses below 1O 6 M 
is very difficult. 

More questionable is the existence of the 
black holes with the masses of hundreds- 
thousands of the solar mass. Determination 
of Mbh in several ULX sources, which are 
hosted by galaxies taken from Messier and 
NGC catalogues, and the discussion of the stel- 
lar mass versus IMBH interpretation can be 
found in Zampieri & Roberts (2009) (see also 
references therein). 

4.3. Non-Active galaxies or Weakly 
Active galaxies 

The best example is of course the Milky Way 
galaxy, with the recent mass determination (see 
Sect. 12. U . The mass in other nearby galaxies is 
known less accurately (for example, triple nu- 
cleus of Andromeda hosts a black hole mass in 
P3 of (1.1 - 2.3) x 10 8 M^. lBender et al. 112005b 
but mass estimate (usually through stellar dis- 
persion me asurement ) is available for hundreds 
of objects dHo 112008b . 

4.4. Active Galactic Nuclei 

The number of known AGN is of order of 10 5 
(in SDSS DR7) and rising due to massive sur- 
veys. The best results are obtained from water 
maser discussed in Sect. 12.11 Generally, mass 
determination for tens/hundreds/thousands of 
objects were perform e d by several authors (e.g 



2002 



Vestergaard & Peterson 
20081 IShenetaL 
20081: iKellyet al. Il2009t IWu 



Woo & UlTV 
2006; Vestergaard et al 
| 2008[lFine et al. 
2009). 

Numerous large surveys supplemented 
with the black hole mass and the Eddington 
ratio determinations (e.g. iHickox et al. 1 2009) 
the available scaling laws from the list dis- 
cussed in Sect. 12.31 The Eddington ratio is 
the frequently determine using another scal- 
ing to derive the bolometric luminosity, for ex- 
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ample from th e optic al monochromatic flux. 
iRichards et all (120061) (see their Fig. 12) anal- 
ysed spectra energy distributions of quasar 
type 1 and gave a formula: 

Lboi = BC A x AL A , (26) 

where the bolometric corrections BC^ = 3.81, 
5.15, 9.26 are for measurements of luminosity, 
AL A , made at A= 1350, 3000 and 5100 A, re- 
spectively. 

It is also possible to use the near-IR 
monochromatic luminosity 

L b0 , = 6.4L(NIR) (27) 

dCaol 120051) or from the broad band X-ray 
lumin osity dHopkins. Richards & Hernquist I 
120071) 

Lboi = (10 - 20)L X (0.5 - 8keV). (28) 

More complex, luminosity- 

dependent corrections can be fou nd i n 
dHopkins. Richards & Hernquist I 2007). 
Hicko xet al. I (12009) used the coefficient in 
Eq. (l2Tb by a factor of two lower than deri ved 
by dHopkins. Richards & Hernquist I I20O7T) in 
order to have consistent results for AGN with 
both IR and X-ray measurements available. 

In heavily obscured Seyfert 2 galaxies or 
radio galaxies the bolometric luminosity can be 
estimated from the [O in] line luminosity: 

L bo i =CL([0 in]), (29) 

where C = 87, 142 and 454 for logarithmic lu- 
minosity ranges of [O ill] line 3 8 - 40, 40 - 42 
and 4 2 - 44, correspondingly dLamstra et al. I 
120091) . 

5. Discussion 

The large number of various methods allow to 
estimate the black hole mass in most astronom- 
ical objects quite reliably, i.e. within a factor of 
a few. Accurate determination of the error of a 
particular measurement is difficult since the er- 
rors are mostly systematic. 

However, there are still classes of objects in 
high need for better mass determination. Those 



are first of all black holes in ULX and in glob- 
ular clusters. The existence of IMBH is an in- 
teresting question both from the point of view 
of the accretion pattern (stellar mass black hole 
hypothesis for ULX sources requires very high 
Eddington ratios and strong beaming). Also 
some other types of active nuclei like weak line 
quasars (WLQ) and radio-loud objects like BL 
Lacs are still in a need for mass determinations 
although some estimates ar e possible even in 
those difficult cases (see e.g. lHryniewicz et aT71 
120091) . 

The mass determination together with the 
estimate of the bolometric luminosity opens a 
possibility of statistical studies in case of mas- 
sive black holes since thousands of objects are 
available. Such a black hole spends most of the 
time in a relatively quiet state. Eddington ra- 
tios cover the range from 10~ 9 to ~ 1, with 
LINERS (and S eyfert 2 galaxies) grouped at ~ 
10~ 4 dHoi|2Q08l), Seyfert 1 galaxies at ~ 10~ 2 , 
and quasars as well as Narrow Line Seyfert 1 
galaxies close to 1. Future evolutionary stud- 
ies should explain the observed ratios of these 
types of galaxies. 
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DISCUSSION 

ANDRZEJ ZDZIARSKI: What is the cur- 
rent range of black hole masses in binaries and 
AGN? 

BOZENA CZERNY: The smallest one cur- 
rently measured is 3.8 solar masses in ob- 
ject J 1650-500 and comes from the paper by 
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Shaposhnikov & Titarchuk. As for the heaviest 
one, they are among the most distant quasars. 
For example, the black hole mass of z — 6.41 
quasar SDSS J l 148+5251 was e stimated to 
be 3 x 10 9 M G dWillot et al. 1120031) . and some 
quasars in Verstergaard (2009) sample have 
black hole masses up to 3 xl0 10 M Q . The black 
hole in OJ 287 mentioned in the context of 
binary black holes, with its mass of ~ 1.8 x 
10 10 M o is also among the largest ones. 

WOLFGANG KUNDT: In Vestergaard et al. 
SDSS plot of the quasars as a function of red- 
shift z, their average masses evolve from a few 
times 1O 9 M , at z <= 5, down to a few 1O 6 M 
at z <= 0.2. How can you downsize this in- 
verted evolution to their predicted growth via 
accretion? 

BOZENA CZERNY: The plot shows only 
very active galaxies. The apparent rise of the 
quasar black hole mass with the redshift in that 
plot, and in the fits by Labita et al. (2009) to 
the quasar mass as oc (1 + z) l M only shows 
that more massive local over-densities evolve 
faster, which is known in cosmology as anti- 
hierarchical evolution. Those large mass over- 
densities becomes quasars at high z, and later 
on pass to a quiet non-active stage, not show- 
ing up in AGN surveys. Smaller over-densities 
evolve slowly, they become active much later 
and never grow up to give large black hole 
masses. For example, Sgr A* is likely to have 
occasional more active stages than seen at 
present but will never become a billion solar 
mass black hole. 



